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We present calculations of the Casimir interaction between a sphere and a plane, using 
a multipolar expansion of the scattering formula. This configuration enables us to study 
the nontrivial dependence of the Casimir force on the geometry, and its correlations with 
the effects of imperfect reflection and temperature. The accuracy of the Proximity Force 
Approximation (PEA) is assessed, and is shown to be affected by imperfect reflexion. Our 
analytical and numerical results at ambient temperature show a rich variety of interplays 
between the effects of curvature, temperature, finite conductivity, and dissipation. 
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Introduction 

Measuring the Casimir force^ has been the aim of an increasing number of experi- 
ments in the past fifteen years (see Ref. 2 and references therein). The comparison 
of these measurements with theoretical predictions from Quantum ElectroDynamics 
have been apphed to put constraints on hypothetical new forces predicted by uni- 
fication models.^'"' Accurate theoretical computations, which account for a realistic 
modeling of experimental conditions, are sorely needed for all comparisons to be 
rehable.^'^ 

The configuration initially considered by Casimir^ is a pair of two infinite par- 
allel plates, perfectly reflecting and at zero temperature, which constitutes a much 
idealized case. The effect of finite conductivity plays an essential role in the accurate 
determination of the force, ^ while the thermal fiuctuations give rise to a remark- 
able interplay with the former effect.* Indeed, in the calculations performed for the 
geometry of two parallel plates, the Casimir force computed from the dissipative 
Drude model turns out to be a factor of 2 smaller than the result obtained from the 
lossless plasma model. 

Moreover, the geometry of two parallel infinite plates is quite particular, for 
example the reflection on the mirrors is specular and the two electromagnetic po- 
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larizations are uncoupled. These features are not present in the general case, and 
the spectral nature of the Casimir effect makes it sensitive to changes in the ge- 
ometry. The nontrivial situation of a sphere and a plane, the configuration of the 
most precise experiments, is depicted in Figure 1: a sphere of radius R is located 
at a distance L from an infinite plane along the ^-direction, L being the distance of 
closest approach and C = L + R the center-to-plate distance. In this configuration, 
not only curvature, but also non-specular refiection, coupling of electromagnetic 
polarizations and finite size are present. 




Fig. 1. Sphere of radius R and a fiat infinite plate at a distance L, separated in the z-direction. 
The center-to-plate distance is C = L + R. 



In the sphere-plane configuration, the Casimir force is usually derived from the 
Casimir energy evaluated in the parallel-plate geometry with an integration over the 
sphere surface. The resulting proximity force approximation (PFA) for the Casimir 
force is^ 

FpFA(i:,i?) = 27ri?^^^^ , (1) 

where is the Casimir energy per unit area in the plane-plane configuration. 
Eq.(l) is expected to provide an accurate description in the limit of small aspect 
ratios ^ (see Refs. 10, 11, 12 for derivations with perfect mirrors at zero tempera- 
ture) . Even if the spheres used in the experiments are much larger than the distance 
from the plane, mastering the error made by this approximation remains necessary 
in order to match the experimental accuracy level. ^'^ Moreover, there is no reason 
why the thermal, finite conductivity, and beyond-PFA corrections should be ex- 
pected to be independent. Therefore, an accurate description of the sphere-plane 
configuration has to take these effects into account simultaneously within a single 
theoretical model. 

In this paper, we present a complete analysis of the Casimir interaction in the 
sphere-plane configuration, using a multipolar expansion of the scattering formula 
for the evaluation of the force. This enables us not only to assess the accuracy of 
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the PFA quantities, but also to investigate the nontrivial dependence of the Casimir 
effect on the geometry, and its interplay with the effects of temperature and optical 
properties of the mirrors. In the first section, we develop the scattering formula for 
the sphere-plane configuration and evoke, in the prospect of numerical evaluations, 
the need for using a finite number of multipoles in the expansion. Section 2 discusses 
the results at zero temperature, first for perfect and then for metallic mirrors. The 
ambient temperature case is finally investigated in Section 3. 

1. Method 

We use the scattering approach^ to compute the Casimir interaction in the sphere- 
plane configuration at ambient temperature. The Casimir free-energy is written as 
a sum over the Matsubara frequencies (n > 0): 

J^ = kBTY,^ndet[X-M{Cn)] , in=n^^^ (2) 

n 

= 7^s(Oe-''^^^'^7^p(Oe-'=«)^ , (3) 

where the primed sum means that the (n = 0)-term is counted for a half. The 
reflection operators of the sphere, 'R-siK)^ ^'^id the plate, TZp{£,), are evaluated with 
reference points at the sphere center and at its projection on the plane, respectively. 
The operator e"'^*^^-''^ accounts for one-way propagation along the z axis between 
these points, separated by the length jC. Thus, the operator A4{^) represents one 
round-trip propagation inside the open cavity formed by the two surfaces. At zero 
temperature, the discrete set of Matsubara frequencies becomes continuous, and the 
sum has to be replaced by a continuous integral over all frequencies in 1R+. 

For a wave of imaginary frequency w = i^, we first introduce the electromag- 
netic multipoles basis | i?,TO,P, ), where £{£ + 1) and m denote the usual angular 
momentum discrete eigenvalues (with £ > 1 and —£< m < £), and P = E,M for 
electric and magnetic multipoles. This set of modes is well adapted to the spherical 
symmetry of TZs. The planar wave basis is however better adapted for the planar 
reflection operator TZp and the two translations operators e'^^^^^, which are diag- 
onal in this representation. We denote those planar modes by the representation 
I k, (p,p), where k is the wave- vector component parallel to the xy-pla.ne, ^ = ± is 
the direction of propagation along the z-axis, and p = TE,TM the polarization. 

We will thus use both sets of electromagnetic modes, by expressing the round- 
trip operator M{^) in the multipole basis, inserting the identity in the planar wave 
basis (see Ref. 14 for a more detailed derivation). The expressions of the non-zero 
diagonal blocks (mi = m2 = m) then read: 




xrp(k,0e-''=^(k,-,p|£2,m,P2> (4) 
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This expression has a simple interpretation when read from right to left: a 
niultipole wave | £2,m,P2) is first decomposed into plane waves (coefficients 
(k, — ,p I ^2,^1,^2)) which propagate towards the plane (factor e~'^^). After re- 
flection by the plane (specular amplitude rp(k, ^)), the plane wave components 
propagate back to the sphere (second factor e~'^^) and are finally scattered into a 
new multipole wave | £1, m, Pi). 

Eq.(4) gives an exact formula for the expression of the Casimir free-energy and 
its various derivatives in the sphere-plane geometry, valid for arbitrary values of all 
the parameters. However, in order to evaluate numerically this quantity, one has to 
truncate the dimension of the operator A4. This can be done by setting a maximum 
value for the quantum number £ < £max- From the localization principle, the value of 
^max required to obtain a numerical evaluation for a given accuracy level is expected 
to scale with the dimensionless parameter ^ = As the frequencies giving the 
main contribution to the Casimir effect scale as ^ < it follows that the required 
^max should scalc as ^. 




Fig. 2. Casimir energy E computed with the scattering formula, normaUzed by the PFA quantity 
-EpFA, as a function of the cut-off ^max, for different values of the aspect ratio 



Figure 2 confirms this effect of the truncation by presenting the Casimir energy at 
zero temperature E, normalized by the PFA result E'pfa- We first consider a sphere 
that is five times larger than the distance to the plane (^ ~ 0.2, circles). When 
the truncation value i^max increases, the result of the numerical evaluation quickly 
converges for the value £max = 20, approximatively. Then, for a twice bigger sphere 
(^ — 0.1, squares), the same convergence is observed, but it is twice slower, such 
that accurate evaluations need i?max ^ 40. This proportional behavior is observed 
as well for even larger spheres (^ = 0.05, triangles). 

From this observation we can conclude that for a given ^max in the numerical 
evaluation, results are only accurate for aspect ratios larger than a minimum 



value, which is inversely proportional to 



This multipolar treatment is then 



February 16, 2012 1:47 WSPC/INSTRUCTION FILE 

benasque_proceedings_canaguier 



The Casimir effect in the sphere-plane geometry 5 

well-suited for situations where the aspect ratio takes intermediate and large values 
!)■ However, the use of great values for fmax in the numerics can also enable 
evaluations in the opposite regime, where ^ takes smaller values. 

2. Results at zero temperature 
2.1. Perfect mirrors 

We start our analysis with the simplest case of perfectly reflecting mirrors at zero 
temperature. Only two length scales are involved in this case, the sphere radius R 
and the distance L between the two objects. The dependance of the Casimir effect 
is then strictly geometrical, and the important parameter is the aspect ratio ^. 

We compare the numerical results for the Casimir energy E obtained from the 
scattering formalism to the energy computed within PFA, .Bpfa- The correction 
factor pe = then a function of ^ only, and does not depend separately on 

L and R. For a vanishingly small value of ^ , the domain of validity for the PFA 
should be recovered and the ratio ps should then go to unity. Assuming a Taylor 
expansion of for small values of ^: 

a usual way to assess the accuracy of the approximated quantities EppA is to esti- 
mate the linear correction coefficient /3e- A similar treatment for the Casimir force 
F = — 1^ or force gradient G = — leads to the linear coefficients Pf and Pg- 
Scalar derivations^'^'^^'^^ give the estimation Pe = \ — ^ — —0.173. An experimen- 
tal study'^"^ has also studied this linear correction term and has put a constraint for 
the force gradient \Pg\ < 0.4. 

In Fig. 3 we present the numerically obtained values for this ratio with red 
squares. They lie below unity, meaning that PFA always overestimates the magni- 
tude of the energy, and indeed tend to unity as the aspect ratio vanishes. In order 
to estimate the linear coefficient Pe from this curve, we use the leftmost points 
(blue circles) to extrapolate the data to unity, with the help of a quadratic fit (solid 
curve). We finally extract the initial slope from this fit, which is the linear correction 
coefficient Pe- 

With this procedure, we obtain Pe — —1.47, a number eight times larger than 
the scalar result, however in good agreement with the results of Refs. 17, 18. As the 
scattering formalism considers the electromagnetic field with the two polarizations 
simultaneously, we conclude that the coupling of polarizations plays an important 
role in the Casimir effect between a sphere and a plane. The corresponding linear 
coefficient for the Casimir force gradient Pg — ^ — —0.49 lies outside the experi- 
mental bound of 0.4 and would seem to be in contradiction with it. However, finite 
conductivity effects play a significant role in the experiment, so we need to consider 
the full theoretical model for real metals, discussed in the next subsection, in order 
to compare with the experimental bound. 
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Fig. 3. Correction factor = for the Casimir energy, as a function of the aspect ratio 

The mirrors are perfectly reflecting at zero temperature. Red squares represent the numerically 
obtained data with the scattering method, the blue circles indicate the points used to compute 
the extrapolating polynomial, presented by a solid line. 



2.2. Metallic mirrors 

The next step in the reahstic description of the experimental configuration is the in- 
clusion of the effects of finite conductivity on the mirrors. This can be done through 
the plasma and Drude dielectric functions, 

expressed at imaginary frequencies uj = z^, where wp is the plasma frequency and 7 
is the relaxation frequency associated to dissipation. The plasma dielectric function 
gpias obtained from Eq. 6 in the lossless limit 7 — 0. 

We then conduct the same procedure as in the previous case but with mirrors 
described by the plasma model^^. We study the Casimir force gradient G in the 
case of a nanosphere with radius R — 100 nm, in order to be able to observe the 
effects of imperfect reflection at short distances {L < Xp = Fig. 4 presents the 
obtained results, for perfectly reflecting mirrors in red and for the plasma model in 
green. 

First, the fact that the two curves do not superimpose means that the accuracy 
of PFA is affected by the introduction of the finite conductivity in the mirrors. 
The two curves do not have the same shape at small distances: while in the case 
of perfect mirrors pc is always a convex function of the curve corresponding 
the plasma model (green curve in Fig. 4) shows an inflection point and becomes 
concave for small values of the aspect ratio. This bending lowers the absolute value 
of the linear correction coefficient, which becomes /3g — —0.2 when using the plasma 
model for the description of the mirrors. This number is back in agreement with 
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Fig. 4. Correction factor pa = 



G 



for tlie Casimir force gradient at zero temperature, as 



a function of the aspect ratio The perfectly reflecting case is presented in red, the case of 
mirrors described by the plasma model in green. Crosses represent the numerical data, accurate 
for > 0.4. The data points used for the fit are specified with circles, and the red dashed-line 
and green solid line are the obtained polynomial fits. The plasma wavelength is Ap = 136 nm. 

the experimental bound |/3gI < 0.4, and is a sign of correlations between the effects 
of geometry and finite conductivity. 

3. Results at ambient temperature 
3.1. Perfect mirrors 

We now include the thermal fluctuations to the vacuum fluctuations of the elec- 
tromagnetic field, through the summation over Matsubara frequencies. In order to 
observe the effect of temperature on the Casimir force, we introduce the thermal 
correction factor^" for the Casimir force -d — -p^, and first consider the case of 
perfect mirrors. A similar factor i?pfa can be derived for PFA quantities. 

The thermal correction factors d obtained for different values of the sphere radius 
R are presented in Fig. 6 with different colors. The first observation is that we always 
have < i?pfAi meaning that the PFA quantities overestimate the contribution of 
thermal photons to the Casimir force. Moreover, a correct treatment of the geometry 
enables us to observe a non-trivial i?-dependence for the thermal correction factor 
•dp- This is is a clear sign of correlations between the effects of temperature and 
geometry, which have been observed for scalar fields. 

For small spheres (i? < 2 /im), we observe that i? < 1 over a large range of 
distances, a feature that never appears for the PFA quantities. In such case the 
thermal photons have a repulsive contribution to the Casimir force, which can be 
associated with the appearance of negative values for the Casimir entropy^^ S = 
~§f ' where J- is the Casimir free-energy. 



February 16, 2012 1:47 WSPC/INSTRUCTION FILE 

benasque _ proceedings _ canaguier 



8 Canaguier- Durand, Lambrecht, Reynaud, Guerout 



1.1 
1.08 
1.06 



t 1.04 

0) 
Q. 

1.02 



0.98 



0.9£ 



R = 0.1 urn 




R = 0.2 urn 




R = 0.5 urn 




R = 1 urn 




R = 2 iim 




R = 5 




R = 10 ym 


t J / I 1 1 III 


PFA 


/ / / 1 1 1 IIt 


— • — analytic 

















L[/xm] 



Fig. 5. Thermal correction factor i9 = -p^y for the Casimir force in the sphere-plane geometry, 
as a function of the distance L. The PFA result i?pfa is recalled by a dashed-line, the solid lines 
are the results computed from the scattering formula, with different colors for different sphere 
radii. The analytical limit of small spheres {R <C L) is drawn as a dotted-line. The temperature is 
T = 300 K. 



This observation is confirmed by the analytical result for the small sphere limit 
{R <§; L) that can be obtained from the corresponding long-distance limit (LD) of 
the Casimir free-energy, introducing the dimensionless parameter v — : 



f 3hcB^ V sinh v + cosh v {v"^ + sinh^ v) 
= -il^'^(-) -th ^(.) = ^-j^3- , (7) 



presented with a dotted-line on Fig. 6. Let us note that the low-temperature limit 
(i? ^ L <C At) of Eq.(7) agrees with the corresponding limit in Ref. 22. 



3.2. Metallic mirrors 

Finally, we include altogether in the model realistic descriptions for the geometry, 
temperature and optical response for metallic mirrors, by computing the Casimir 
force F for metallic mirrors at ambient temperature. We concentrate on the effect 
of dissipation on the Casimir force, by studying the ratio pui^d of the force for 
dissipative and non-dissipative metals. In the case of two parallel plates at non-zero 
temperature, as well as for the derived PFA quantities, this ratio shifts from unity 
at small distances to the value 2 at large distances. 

In Fig. 6 we present this ratio as a function of the distance L for different sizes of 
the sphere. The curves differ from the one obtained from PFA quantities, which is 
a sign of correlations between the effects of geometry, temperature and dissipation. 
Moreover, the long-distance limit of this ratio is no more 2 as in the plate-plate 
configuration, but a rather an increasing function of the dimensionless parameter 
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Fig. 6. Ratio of the Casimir force for lossless plasma and dissipative Drude models with 

respect to the distance L. The PFA result is recalled with a dashed-line, the solid lines are the 
results computed from the scattering formula, with different colors for different sphere radii. The 
temperature is T = 300 K, the plasma wavelength is \p = 136 nm and the length scale associated 
to the dissipation is = 22L£ = 34 ^m. 



a = 27rj^, that can be derived analytically:^ 

, 3 / 1 cotha\ 

/(«)=o 1 + — ' 8 

2, \ J 

taking values between 1 and |. This observation hints at the existence of a strong 
interplay between the geometric, thermal, and dissipative corrections. It also means 
that the gap between the results obtained from plasma and Drude model is much 
smaller than expected from the PFA quantities. 



Conclusions 

In this paper we have investigated the Casimir interaction in the sphere-plane con- 
figuration, which can be seen as a simple example for studying the dependence of the 
Casimir effect on the geometry. The scattering approach allows to treat exactly this 
configuration, taking into account both the non-zero temperature and the metallic 
nature of reflectors. 

First, we have observed that the coupling of polarizations plays a role in the 
Casimir effect, a complete electromagnetic treatment is thus necessary to get reliable 
evaluations in this geometry. The inclusion of a realistic description for the material 
properties also affects the accuracy of the quantities obtained from the PFA, and 
it results in a much smaller linear correction coefficient j3Q which lies within the 
experimental bound. 

At ambient temperature, the correlations that we observe between the effect of 
geometry and temperature can lead to a repulsive contribution of thermal photons 
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to the Casimir force, and to negative values for the entropy. Finally we discuss 
the effect of dissipation on the materials, and the long-distance ratio between the 
force obtained from plasma model and Drude model, which is 2 in the parallel-plate 
geometry, is reduced to, at most, 3/2. 
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